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Since  final  Report  BAJA37-69-C-1223,  the  following  papers 
have  been  published  on  work  supported  partially  by  this  Contract: 

(i)  Interaction  between  point  defects  in  metals,  by  N.  H.  March 
and  J.  S«  Rousseau,  Crystal  Lattice  Defects,  2,  1-1:6,  1971 

(ii)  Partition  function  of  a  disordered  system,  by  J.  S.  Rousseau, 
J.  C.  Stoddart  and  N  H.  March,  J.  Fhys.  C.,  U,  L59-L&1,  1971 

(iii)  The  Dirac  and  t  matrices  for  spherically  syBtnetric  potentials, 
by  J.  S.  Rousseau,  J.  Pfcys.  C.,  It,  L391-L39U,  1971 

Furthermore,  an  invited  paper  was  given  at  the  Battelle 
Colloquium,  Seattle  last  Summer  on  1  Computer  Simulation  of  Lattice 
Defects*  and  will  be  published  in  the  Proceedings  of  that  meeting. 
This  is  entitled 

•Ranzaki  forces  and  electron  theory  of  relaxation  round  defects* 
by  N.  H.  March  and  J.  S.  Rousseau. 

Ihe  substance  of  this  latter  paper  is  reported  in  Part  a  of  the 
present  report. 
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Part  B  continues  our  work,  begun  by  Rousseau,  Stoddart  and 
March  (1970)  and  continued  in  (ii)  above,  on  electron  states  in 
disordered  systems,  by  considering  specifically  the  electrical 
resistivity  of  liquid  metals  in  the  strong,  as  well  as  the  weak 
scattering  regime.  While,  for  weak  scattering,  wa  regain  the 
Ziman  formula,  we  expect  our  theory  also  to  apply  to  strong  scatterers, 
HJra  the  liquid  alkaline  earth' metals.-  - '  Of  particular  intorewt 
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here  are  the  metals  calcium  and  barium,  for  which  electrical 
transport  measurements  have  become  available  very  recently 
(Van  Zytveld,  Enderby  and  Collitigs,  1972).  These  workers 
conclude  that  no  simple  modification  of  the  pseudopotential  theory 
can  exp?tain  their  experiments.  We  report  in  Part  B  therefore 
a  basic  theory  which  is  capable  of  dealing  with  strong  scattering 
such  as  we  encounter  in  these  liquid  metals.  Calculations  of 
the  single-centre  scattering  for  calcium  are  in  progress,  in 
order  to  br;.ng  our  theory  based  on  the  inverse  transport  coefficients 
into  contact  with  experiment. 


PART  A 


1.  Intooduction 


Especially  in  metal3,  there  has  been  a  good  deal  ol  interest 
in  the  theory  of  wiarged  defects;  for  example  a  vacancy  in  Qu, 
or  Mg  a  tens  in  dilute  concentrations  in  Id.  So  far,  however,  no 
very  careful  account  of  the  detailed  relaxations  of  the  lattice 
round  the  defect  or  impurity  site  has  been  incorporated  in  such 
electron  theory  calculations,  though  attempts  have  been  made  to 
include  some  account  of  relaxations  by  modifying  the  Friedel  sua 
rule,  which  is  frequently  used  in  defect  calculations  to  simulate 
self-consistency. 

Thu  electron  theory  of  defects  in  metals  has  proved  to 
make  interesting  predictions  about,  for  example,  charged  defect 
interactions,  especially  in  polyvalent  metals,  which  appear  to 
agree  in  general  terns  with  experiment  (see,  for  example,  the 
review  by  two  of  the  present  writers  (March  and  Rousseau  1970) ) . 

But  a  substantial  body  of  work  also  exists  in  which  a 
markedly  different  philosophy  is  adopted.  Here,  relaxation 
effects  are  studied  as  tne  prime  objective,  usually  on  the  basis 
of  pair  forces.  Unfortunately,  such  an  approach  often  fails 
to  include  a  proper  account  of  electron  redistribution  caused  by 
the  introduction  of  defects  or  impurities  into  the  lattice. 

Though  this  is  a  serious  limitation,  nevertheless  the  approach  has 
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p^ctical  merit  in  that  it  allows  local  relaxations  round  defects 
to  be  estimated,  in  a  way  which  matches  these  on  to  the  correct, 
long-range  lattice  displacements  given  by  elasticity  theory,  at 
large  distances  from  the  defect.  Indeed,  quite  a  body  of  work 
ic  now  available  giving  us  estimates  of  the  displaced  positions 
of  lattice  atoms  about  a  vacancy,  or  a  divacancy,  in  a  number  of 
metals. 


Such  calculations,  mace  with  plausible  pair  forces,  show 
that  even  when  a  vacancy  in  created  in  open  body-centred  cubic 
structure  like  Na  metal,  the  atonic  movements  are  really  quite 
small.  This  is  a  fortunate  circumstance  and  had  been  anticipated 
by  Kanzaki  (1957;  sea  also  Matsubara,  1952)  in  his  method  of 
lattice  statics.  Here-  the  idea  is  to  simulate  the  effect  of 
a  vacancy  or  impurity  in  a  lattice  by  regarding  the  originally 
perfect  lattice  as  strained  until  it  takes  up  the  relaxed 
configuration  appropriate  to  the  equilibrium  atomic  positions 
in  the  defect  lattics.  Obviously,  external  forcos  must  be 
applied  to  hold  the  other  atoms  in  tfceir  displaced  positions. 

We  shall  refer  to  these  external  forces  as  Kanzaki  forces  and 
to  the  3trainod,  but  otherwise  perfect,  lattice  as  the  Kanzaki 
lattice.  Kanzaki  assumed  that  all  the  displacements  u.  , 
say,  from  the  perfect  lattice  positions  _£ .  were  sc  small  that 


one  need  work  only  to  first-order  in  u  . 

V 


Thus,  on  the  one  hand,  ve  have  the  electron  theory,  in 


jk 

3 


! 


I" 

-3- 

which  a  basic  quantity  characterizing  the  defect  is  the  localized 
charge  it  displaces  (as  a  specific  example  to  be  referred  to  again 
later,  see  the  calculations  of  Stott,  Baranovsky  and  March,  1970, 
on  the  charge  displaced  round  vacancies  in  close-packed  polyvalent 
metals)  and,  on  the  other,  the  method  of  lattice  statics,  with  the 
strained  Kanzaki  lattice  held  in  the  relaxed  final  configuration 
by  external  forces. 

Though  the  present  report  is  concerned  with  the  electronic 
structure  of  defects,  the  idea  of  the  Kanzaki  lattice  will  play 
a  central  role  in  setting  up  the  theory.  Specifically,  this 
allows  ua  to  develop  a  means  of  treating  the  difference  between 
the  charge  displaced  by  a  given  defect  potential  in  the  unrelaxed 
lattice  and  that  displaced  by  the  same  defect  potential  inserted 
in  toe  Kanzaki  lattice.  Sane  p.'eliminary  results  for  a  vacancy 
in  Cu  metal  are  referred  to,  by  way  of  illustrating  the  approach. 

It  can  be  objected  that,  to  set  up  the  Kanzaki  lattice, 
we  must  know  toe  relaxations  at  the  outset.  This  is  true,  and 
we  must  use  some  starting  estimates  such  as  those  referred  to 
above  from  pair  potential  studies.  However,  once  the  charge 
displaced  in  the  Kanzaki  lattice  by  the  defect  potential  is  known, 
we  are  in  a  position  to  calculate  the  electric  field  at  toe 
ionic  positions  in  toe  Kanzaki  lattice,  and  hence  from  the 
Feynman1 s  theorem,  the  Kanzaki  forces.  If  these  ai.^e  not 
consistent  with  the  original  displacements  ,  Kansaki’s  method 
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can  then  be  used  to  generate  a  new  set  of  displacements.  The 
new  displaced  charge  must  then  be  found.  dearly,  we  can,  in 
principle,  iterate  until  a  self-consistent  result  is  obtained c 


Though,  therefore,  at  least  in  principle,  the  theory 
given  here  can  solve  the  problem  of  the  electron  distribution  round 
a  defect  in  a  correctly  relaxed  lattice,  there  are  a  number  of 
practical  issues  to  be  resolved,  which  will,  almost  certainly, 
involve  major  electronic  computation.  Of  these,  a  prime 
problem  is  that,  to  obtain  the  electronic  structure  of  a  defect 
centre,  wa  must  have  a  great  deal  of  information  about  the  perfect 
lattice.  We  do  not,  as  yet,  have  that  information  in  a  suitably 
explicit  form  for  any  metal  but,  as  we  shall  show  in  section  3 
below,  a  good  deal  of  it,  though  not  all,  can  be  got  frcm  well 
established  methods  in  the  band  theory  of  perfect  crystals.  The 
second  difficulty  resides  in  the  fact  that  methods  of  solving 
the  problem  of  scattering  off  a  defect  potential  in  the 
unrelaxed  lattice  are  still  troublesome  to  apply.  However,  when 
both  the  perfect  crystal,  and  the  defect  potential,  can  bo 
constructed  from  muffin  tin  potentials,  we  have  the  rigourous 
one-electron  theory  of  Beeby  (1967),  while  for  three  classes  of  non- 
muffin  tin  defect  potentials  we  have  available  the  explicit 
approxL  ntiotfe given  by  Sfcoddart,  &rch  and  Stott  (1969).  A 
third  difficulty  is  that  we  do  not  know  precisely  what  potential 
to  take  to  describe  the  defoct.  But  at  least,  we  can  give 
an  operational  definition  of  what  that  potential  must  be  and 
wo  do  this  in  section  h. 
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We  concise  this  part  of  the  report  with  a  brief  discussion 
of  the  way  defect  energetics  might  be  studied  within  the  present 
framework. 
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A  convenient  starting  point  from  which  to  explain  the  basic 
ingredients  needed  in  the  electron  theory  is  to  study  how  the 
electron  density  p  Q(r)  in  the  perfect  crystal,  which  is  accessible 
•■ to  experiment  via  the  intensities  of  X-ray  scattering  at  the  Bragg 
reflections,  is  changed  when  we  strain  the  crystal  to  fora  the 
Kanzaki  lattice,  and  we  turn  immediately  to  the  solution  of  this 
problem. 

2.  Electron  density  in  Kanzaki  lattice 


We  start  with  a  perfect  crystal,  with  atoms  at  lattice  points 
We  simply  strain  the  lattice,  introducing  neither  foreign  atoms,  nor 
vacancies. 


The  method  we  use  to  generate  the  electron  density  in  this 
Kanzaki  lattice  is  that  used  by  Jones  and  March  (1970,  referred 
to  subsequently  as  JM)  in  their  theory  of  lattice  dynamics.  Ute 
theory  is  immediately  useful  within  the  framework  of  the  Kanzaki 
method  of  lattice  statics,  provided  only  that  the  displacements 
from  the  sites  ^  ,  say  u_^  ,  are  sufficiently  small  so  that  we 
need  work  only  to  first-order  in  the  displacements.  This  is  a 
basic  assumption  of  our  approach,  and  it  will  have,  of  course,  to 
be  checked  for  consistency  in  any  application.  As  >re  stressed  above, 
however,  even  with  a  vacancy  in  an  open  structure,  it  is  found  in 
practice  that  relaxations  are  often  quite  a  small  fraction  of  the 
lattice  parameter,  and  the  present  approach  should  be  widely  applicable. 
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Iet  the  change  in  electron  density  when  the  lattice  is 
•strained*  bop^(r).  Then,  to  first-order  in  the  displacements,  the 
theory  of  JM  allows  us  to  write 


p  (r)  =  2  u,  .  B(r  -  £). 
1  e  *  “ 


(2.1) 


The  vector  R  satisfies  t n  integral  equation  given  by  JM  which 
involves  knowledge  of  the  exchange  and  correlation  energy  ox'  the 
inhomogeneous  gas  of  electrons  in  the  crystal  (see  Appendix  .  for 
some  relevant  details.)  However,  it  is  important  to  note  that  R(r) 
evidently  also  deteisdnos  tha  gradient  of  the  perfectly  periodic 
lattice  densityp  (r)  through 


?Po  (r)  =  2  R(r  -  £), 

£  ~ 


(2.2) 


Information  about  R(r)  in  Fourier  transform  is  therefore 
available  at  the  lattice  vector*  K.  Here,  in  fact  is 
determined  uniquely  by  the  intensity  of  X-ray  scattering  at  the 
Bragg  reflections,  as  discussed  in  JM.  In  particular 


%  *  1  £  <k 

where  is  the  Fourier  component  at  K  of  the  charge  density 


(2.3) 


Thus,  it  is  dear  that  if  we  know,  from  say  a  pair  potential 
study,  a  first  approximation  to  the  u^’s  and  if,  as  discussed  above, 
these  relaxations  turn  out  to  be  sufficiently  small,  we  can  generate 
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the  density  PQ(r)  +  P1 (r)  at  any  point  in  the  Kanzaki  lattice 
from  a  knowledge  of  the  quantity  R  characteristic  >f  the  perfect 
lattice. 

As  an  example  to  illustrate  this  method,  we  shall  consider 

in  Appendix  3  sane  preliminary  numerical  results  for  R  in  Cu  metal, 

from  which  P  ^ (r)  for  a  vacancy  in  Cu  metal  can  be  estimated  using 

the  displacements  u  calculated  by  Tewordt  (1958). 

“■  £ 

However,  it  is  clear  that,  though  p  ^ (r)  is  a  basic  quantity 
reeded  in  the  theory,  a  gooa  deal  more  information  is  required  in 
order  to  generalize  the  approaches  of  Beeby  (1967)  or  Stoddart, 

March  and  Stott  (1969)  to  deal  with  perturbations  in  the  Kanzaki 
lattice.  We  shall  therefore  turn  next  to  discuss  the  Green  function 
or  density  matrix  describing  the  Kanzaki  lattice.  The  discussion 
of  the  density  above  is,  of  courss,  a  special  case  of  this  more 
general  treatment.  We  want  to  emphasize  though  that  the  above 
discussion  is  exact  to  first-order  in  the  displacements  if  the 
properties  of  the  perfect  lattice,  as  summarized  in  R(r),  are 
known. 
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3*  Green  function  and/or  density  matrix  in  Kanzaki  lattice 

We  can,  as  ve  shew  below,  set  up  the  Green  function  of  the 
Kanzaki  lattice,  from  that  of  the  perfect  lattice,  GQ(r  r‘E)  say, 
the  latter  being  constructed  such  that  the  exact  electron  density 
p  o  (r  r  Ef),  with  Ef  the  FemL  energy,  is  correctly  incorporated. 

Below  ve  shall  work  with  the  canonical  density  matrix 
CQ  (r  r'P)  for  the  perfect  lattice,  defined  by 

C  (r  r’B)  =  Vto*  (v)  ilr  ,  (r‘)  “3  EY(k)  , 

° -  v-'  Yk  ;  e  Y  (3.1) 

ifk 

There  the  HLcch  wave  functions  (r)  and  corresponding  eigenvalues 
Ey  (k)  are  generated  by  the  periodic  potential  V^(r)  vhich  gives 
by  definition  the  exact  ground-state  density  in  the  perfect 
unstrained  lattice.  If  ve  take  the  Laplace  transform  of  CQ 
with  respect  to  0,  and  call  the  transformed  variable  -E,  then 
ve  regain  the  Green  function  0Q  (r  r*  E)  immediately, . so  that 
vhether  ve  work  vtith  C  or  G  is  purely  a  matter  of  mathematical 
convenience. 


Following  JM  ve  can  vrite  for  the  strained  lattice,  a 


perturbing  potential  A  V(r)  having  the  form 
A  V(r)  =  V  u.  .  p  (r  - 


>)  =  ^  He  *  E  (£  -  £) 


(3.2) 


Trvn  JH, (see  also  Appendix  1 )  «e  know  that 
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R(r)  =  yp(xf)  F(r  r1)^’  (3.3) 

whore  F  is  a  one -body  response  function  given  lay  (cf  Appendix  1) 

--  =.•  2ReTG.(r  r'E)  Bo,*  (r  r'E)  i  (3.U) 

3E  L°  3E  J 

where  P  (r  r‘E)  is  the  Dirac  density  matrix  of  the  perfect  lattice. 
Once  the  periodic  potential  Vp  is  known,  P  Q  and  Gq  can  be  calculated 
purely  fraa  the  Bloch  functions  %  (r)  and  the  energies  E  Yk 

It  is  clear  fraa  (3.3  )  that  knowledge  of  the  one-bc  ▼ 

response  function  F(r  r»)  plus  P(r)  is  equivalent  to  knowing  R(r). 

To  generate  the  density  matrix  in  the  Kanzaki  lattice,  we  now 
need  to  L'olve  the  Bloch  equation 

HC  =  -  9C  (3.5) 

.  ap 

where 

H  =  -g  V2  +  Vp(r)  +  AV(r)=  HQ  +  AV,  (3#6) 

subject  to  the  usual  delta  function  boundary  condition  C(r  £c  q) 

m  -  t  )  expressing  the  completeness  of  the  eigenfunctions. 

— o 

But  Hq  generates  the  density  matrix  0Q  and  A  V  is  a  pert’irbation, 
frcm  oqn.  (  3.2  ),  provided  the  u^  s  are  small.  Hius,  we  can 
write,  to  first-orderin  AV  (cf  March,  Young  and  Sampanthar,  1967) 


:.  )^'  -«<«^  *  -^f*V  *T  ~  ^  *'*  ?  • 


li.  Solution  for  defect  lattice 


We  now  introduce  the  defect,  and  to  be  specific,  we  shall 
assume  ws  create  a  vacancy  at  the  origin.  Let  us  suppose  that 
the  electron  density  p^(r)  *  p  (r)  +  p  ^  (r)  in  the  Kanzaki 

lattice  changes  to  p^(r),  the  density  in  the  final  state  containing 
the  relaxed  vacancy. 


Ij.l  Operational  definition  of  defect  potential 


We  tacitly  assumed  above  that  we  could  generate  the  exact 
ground-state  density  PQ(r),  as  observed  say  in  X-ray  scattering, 
from  a  one-body  periodic  potential  Vp(r):  i.e. 


P0(e)  = 


:  l_  <4  (E)  %  (i)  • 

(occupied) 


(i»- 1) 


This  essentially  fellows  from  the  considerations  of  Kohn  and 
Hohenberg  (196U)  and  Kohn  and  Sham  (1565).  An  operational 
procedure  to  construct  V  (r)  from  a  given  P  (r)  has  been 
discussed  rather  fully  by  Beattie,  Stoddart  and  March  (1971). 


The  same  argument,  in  essence,  now  enables  us  to  define 
operationally  the  defect  potential.  Wo  wish  to  find  a  potential 
Vj(r)  which,  when  added  to  the  one-electron  Hamiltonian  H  in 
eqn.  (  3.6  ),  yields  the  exact  final  state  density  Pf(r). 
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We  want  to  stress  that  V^(r)  is  a  potential  to  be  added 
to  H,  not  to  Hq.  Thus,  V^(r)  is  generated  in  a  system  in  which 
the  ionic  configurations  differ  only  by  the  removal  of  ths  ion 
at  the  origin.  If  we  had  defined  the  defect  potential  as  a 
charge  fieri  flQ,  all  the  ions  would  have  soved  and  Vj(r)  would  bo 
a  much  more  complex  object. 

For  the  vacancy,  especially  in  polyvalent  metals,  we  must 
not  expect  that  the  effect  of  7^(r)  can  be  treated  as  a  perturbation. 
The  Bloch  equation  (3.5  )  with  H  replaced  by 

Hf  -  Ho  +  AV  +  Vd  (r) 
can  again  bo  written  as  an  integral  equation 

Gf  (r  r»  p )  «  C^(r  r«  p  ) 

cipl  jar"  CR  (r  r,,p-P1)VdCrw)  Cf(r,,r'P1) 

whore  now  the  object  which  is  required  appears  also  on  the 
right-hand  side  of  eqn.  (  h.j  ).  In  principle,  knowing 
frem  eqn.  (  3.7  ),  we  can  obtain  from  eqn.  (  ^,3  ),  with 

an  assumed  defect  potential.  In  practice,  an  iterative  scheme 
(cf .  Hilton,  March  and  Curtis,  196  7  )  would  have  to  be  used 
and  the  procedure  is  certainly  very  lengthy  and  somewhat  troublesome. 
We  shall  therefore  consider  below  two  approximate  methods  which 
allow  us  to  solve  (or  in  the  second  method  to  circumvent)  the 
integral  equation  '  U.3  )« 


0».2) 

CUo) 


«WSS  V’ ^CT<s^  ’’ 


Hie  idea  has  already  been  made  plain.  Physically,  we 
anticipate  that  the  diagonal  difference  Cf(r  r  p  )  -  0^(t  £  p  ) 


H»  can  be  usefully  split  into  two  parts,  aCuq  -  Cu(r  r  p  ) 


-  C  (r  r  0  ),  there  C  is  the  HLoch  density  obtained  by  introducing 
the  defect  potential’  into  the  perfect  lattice  with  HLoch  density  CQ, 
plus  another  term  6  C  taking  relaxation  into  account. 


U.2  Correcticn  to  displaced  charge  as  calculated  in  unrelaxed  lattice 


A  rough  approximation  to  estimate  &C  can  be  given  as  follows. 


Write,  with  the  assumption  that  V.Qr)  varies  slowly  in  space 
(cf.  eqn  (3.1  )  ) 


Cf(r  r  0)  =  (r  r  p)  e~Va(^ 


(U.U) 


Cu(£  £  3)  =  C0(r  r  0)  e 


-PV-(r) 


(h.$) 


Thus  we  have  from  the  definitions  given  above 


icnc 


(li.o) 


■  "{•■"*  -  } 


(1*.7) 


Therefore  it  follows  that 


V 
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6C  =  e~^d  -ij  - 


■  -1]  “KO. 


(h.  8) 


Eero  then,-  wo  have  a  rough,  hut  quite  practicable,  way  of 
estimating  tha  correction  to  the  displaced  diarge  in  the  unrelaxed 
lattice,  due  to  relaxation,  since  A  is  given  explicitly 
by  eqn.  (  3.7  ). 


In  terms  of  densities,  we  could  i^teraatively  write 


Pf (r  E)  j  pK(r,  E  -  Vd(r)  ) 


(U.9) 


and 


p„(s  E>  =  Po  <£-  E  '  V-2)> 


•  rO 


(ii.10) 


We  stress  that  we  only  use  these  forms  to  estimate  the  correction 
to  the  displaced  charge  APU0  =  Pu(£  E)  -  P0(£  E)  due  to 
relaxation.  Then  we  have 


APfX  “  Apuo  =  ®P 

=  PK(r,  E  -  Vd(r)  )  -  pK(r  E) 

~  Pc(£>  E  ”  ^d(r)  )  +  pQ(r  E) 

AV  ( v*  ^Pf?  fr>  t 


(U.li) 


=  Jd^  AV  (r^ )  P  (r  E  -  7d(r)  )  -F(r  E)  j 


which  follows  from  the  definition  of  the  response  function  F.  Hare 


then  is  a  first  approximation  to  correct  the  displaced  charge  in 

* 

the  un relaxed  lattice.  Since  &  )  is  given  by  eqn.  (  3.2  ), 

it  is  elear  that  6p  is  0  (u  ),  as  required. 

■“*  ■o 

We  want  at  this  stage  to  cement  further  on  the  result 
(U.ll  )•  Through  the  potential  A  V(r),  determined  by  perfect 
lattice  properties  plus  the  u,  »s  there  is  an  obvious  link  with 
lattice  dynamical  properties.  The  displaced  charge  Ap  uo 
due  to  Vj(r)  inserted  into  the  unrelated  lattice  i3  to  be 
corrected  by  the  appropriate  fora  (  U.H) :  this  reduces  to 
sero,  as  it  must,  when  Vj(?)  ■  0. 

Around  the  defect,  however,  wo  need  to  know  F(r  E) 
as  a  function  of  E,  in  order  to  evaluate  eqn.  ( U.ll  ) .  As 
can  be  seen  from  eqn  (  3»U  ),  this  is  a  problem  in  band  theory 

which,  for  a  given  periodic  potential  V^Cr),  is  soluble  by 
existing  methods,  though  it  will  be  an  extensive  computational 

o 

task  Information  about  the  response  function  F(r  r^  E) 
will  lead  to  progress  in  the  field  of  lattice  defects;  another 
bonus  one  could  gain  then  would  be  to  map  out  the  anisotropy 
of  the  displaced  charge  due  to  a  weak  perturbation,  such  ao 
that  due  to  Be  in  Id. . 


A  more  precise,  but  far  more  complicated,  way  of  correcting  the 
displaced  charge  is  given  in  Appends  1 . 


\ 


!, 


I- 


-18- 


t  is  the  t  matrix  for  a  single  muffin  tin  potential,  whils  all 
the  structure  dependence  is  in  G(m)for  which  Beeby  gives  explicit 
expressions.  The  integration  in  eqn.  (li-lf>  )  is  over  the 
Brillouin  zone  of  volume  t  ,  while  R^(r)  is  as  usual  the  radial 
wave  function.  All  the  information  needed  in  eqn.  (  li- 15>  ) 
is  accessible  through  a  KKR  band  structure  calculation. 


By  comparing  eqn.  (  Ii.lli  ),  used  at  the  Fermi  energy 
E  *  F.f,  wJ.th  the  semi-empirical  procedure  based  on  eqn.  (  2.1  ), 
it  will  be  interesting  to  aee  whether  the  approximation  (U.lli  ) 
is  sufficiently  accurate  to  allow  A  V(r)  to  be  found  explicitly. 

If  so,  we  have,  of  course,  through  eqn.  (ii.lli  )  the  energy 
dependence  of  p  ^  (r  E)  required  in  the  evaluation  of  eqn.  (U.U  ), 
which  in  t-ie  approximation  equivalent  to  eqn.  (  Ii.lii  )  becomes 

6p  =  AV(r)  jjj  (r  E^,  -  V^Cr))-  o(r  Ef)  ..  (U.l6) 

Unfortunately,  though  numerical  calculations  are  in  progress, 
results  are  not  as  yet  available  for  p . 

later,  when  we  have  more  knowledge  of  the  basic  response 
functions,  refinement  may  be  carried  out  via  the  method  of 
Appendix  2,  should  it  prove  necessary. 


lb 


W 


If  we  can  make  the  approximation  that  A  V(r)  i3  slowly 
varying  then  we  find  eqn.  (  U-lii  )  and  if  we  apply  this  at  the 
Feted,  level,  we  have  evidently  the  approximate  result 


Pi.  (H  E)  ~  Pl(r  Ef)  0Q(r  B) 

oo(£.  Bf ) 


(1.17) 


Thus,  in  this  approximation  we  simply  scale  the  displaced  charge 
below  the  Fermi  level  E^.  in  the  Kansaki  lattice  with  the  local 
density  of  states  in  tne  perfect  crystal. 


Using  the  first  approximation  for  6p  that 
aPkE)  =  Pl(r  E  -  Vd  (r))  -  p1:(r  E) 


(fc.J*) 


we  see  that  this  takes  the  explicit  form 

6p(r  E)  -  pi  (r  E)  j"o0(r  E  -  Vfl) 

°o'fe  V 

and  hence,  at  the  Tenai  level 


>o<£  s)l 
i0(r\r  (u*] 


6p(il)  =  P^rEf)  [o0(£  Ef  -  Va)  -  1  J  (h  20) 

%  (S.Ef> 

Thus,  given  the  local  density  of  states  oQ  (r  E)  in  the 
perfect  crystal,  from  the  band  theory  calculation,  th«-«  defect 
potential  V^(r)  and  the  displaced  charge  in  the  Kanzakl  lattice 
at  the  Fermi  level,  we  can  calculate  the  correction  to  the  charge 
displaced  by  in  the  unperturbed  lattice.  Of  course,  we  have 
only  circumvented  the  many-electron  nature  of  the  problem  by  using 
Oqn*  (2.1  )  with  R  (r)  found  empirically  frort  the  Bragg  reflection 

intensities.  We  stress  that  wo  can  include  both  aany-electron 
effects  and  many-bcd y  forces  via  the  presence  of  p  ^  (r)  in  the 
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result  (  ii»20  ).  But  of  course  we  are  still  involved  in  the 
assumption  that  V^  is  slowly  varying  and  wo  expect  that  refinement 
of  the  theory  will  lead  to  a  form. 


6p  (r) 


(U-21) 


vjr*~j+  \ 


Thi  square  bracket  in  eqn  (  ^.20  )  is  a  first  approximation  to 
the  modified  function  f. 


U.3  Local  density  of  states  for  Mg  in  Id. 

Dr.  J.  E.  Inglesfield  (private  communication)  has  recently 
made  prelirainary  calculations  of  the  local  density  of  states  o(rE) 
for  Mg  in  Li  and  v-3  wish  to  conclude  this  section  by  mentioning 
these  results  briefly. 

The  results  are  shown  for  two  energies  \  Ej  and  E^  in  ELgs 

1  and  2.  The  variable  used  is  r/r  ,  where  r  is  the  radius  of 

s  s 

the  LL  atomic  cell  (3-265  a.u. ).  The  Ibrad.  energy  of  la  (0.173 
a.u. )  has  boon  used. 

The  dashed  curves  show  the  original  local  densities  of  states 
for  the  pure  materials,  at  each  energy,  while  the  solid  curves  show 
tho  local  density  of  states  for  a  Mg  impurity  in  Id. 

As  Inglesfield  has  emphasized,  something  like  von  lane’s 
theorem  is  being  recovered  as  we  approach  the  Fermi  energy. 


It  is  quite  clear  that  ..j  can  use  the  local  density  of 
states  for  Li  thus  obtained  to  calculate  the  function  in  the 
square  brackets  in  eqn. (  h>20  ).  Calculations  are  currently 
in  progress,  in  conjunction  with  Drs.  Bullough  and  Perrin,  to 
obtain  o  o(n  E)  firm  Beeby‘s  method  for  Ii,  as  well  as  the 
displaced  charge  round  an  unrelaxed  vacancy  in  this  method. 
Eqn.  (  h'20  )  can  then  be  used  to  correct  for  relaxation. 


; 

i 


5*  ■&&*$&*&** 3?  rj^AtfS'k'ap^ «?Jc~  *~i  •&*%£?%&&&&>  T*-V£j 
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5*  Electric  field  in  Kanzaki  lattice 


Let  us  now  consider  how  we  can  calculate  the  forces  acting 
on  the  ions  in  the  Kanzald.  lattice  from  the  electron  theory  developed 
above.  Again,  to  be  specific,  we  will  consider  the  monovacancy, 
in  which  case  a  single  atom  has  been  removed  from  the  origin.  Clearly 
in  this  defect  lattice,  there  are  no  forces  acting  on  the  atoms, 
since  the  relaxed  configuration  is  in  equilibrium.  This  implies, 
from  Feynman’s  theorem,  that  the  total  electric  field  E  (r)  in 
the  relaxed  defect  lattice  is.  identically  zero  when  r  +  u  «  . 

*  V 

Otherwise,  there  would  be  forces  acting  on  the  nuclei.  This 
electric  field  is  created  hy  an  electron  donsitypo(r)  +  p  (r) 

+  Pd(r),  plus  the  fields  of  the  nuclei. 


Now  ve  put  back  the  atom  at  the  origin,  in  our  vacancy 
example.  The  electron  density  in  this  Kanzaki  lattice  is 
P  Q(r)  +P1  (r).  The  electric  field  l^GZa^  (r)  *  e  (r) 
in  this  strained  lactice  is  evidently  determined  by  the  electron 
density  f^(r)  +  (r)  plus  the  nuclear  configuration.  It  is 

then  dear  that  the  resultant  ele etric  field  acting  at  the 
nuclei  must  be  equivalent  to  that  duo  to  the  difference  between 
the  electron  densities  in  the  Kanzaki  lattice  and  in  the  relaxed 
defect  lattice  and  to  the  difference  between  the  nuclea- 
configurations. 


Immediately,  for  the  nonovacancy,  we  see  that  this  field 
must  be  determined  by  the  electrostatic  potential  due  to  the 
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displaced  electron  charge  round  the  vacancy,  plus  that  due  to 
the  absence  of  the  ion  at  the  origin.  This  electron  density  is 
given  by 

pf  -  pK  =  pate) 
f  Apuo  +  [pi  <£»  Ef  ■  vd(^J  “  Pi(£  Vj  (*‘ 


i) 


and  it  is  this  density,  together  vith  the  field  of  the  Absent* 
ion,  ^hich  determines  the  electric  field  at  the  ions  in  the 
Kanzaki  lattice,  and  hence  the  Kanzaki  forces.  Fran  our  knowledge 
of  potentials  and  displaced  charges  in  free  electron  metals,  it 
seems  that,  at  large  r,  by  Taylor  expanding  eqn.  (U-16  )  in  toms 
of  Vd,  ftp  will  be  both  smaller,  involving  the  product  of  A  V$r) 
and  Vd>  and  shorter  range  than  Apuo  ,  though  this  will  have  to 
be  verified  by  detailed  numerical  claculatione.  Thus  at  large  r, 
we  rhall  tentatively  assume  ^Puo  dominate  in  eqn.  (  5.1  ). 


5-1  Sane  preliminary  results  from  free  electron  model 

However,  we  can  readily  estimate  the  field  due  to  this  term 
Apuo  and  it  is  worth  recording  the  result  here  by  way  of 
illustration.  We  have,  from  the  Poisson  equation 


V20  =~4*  Apuo  (5*2) 

where  <f>  is  the  electrostatic  potential.  In  the  case  when 
Apuo  is  spherically  syaxaetrlc,  we  find 

r2  d£  =  [r  4*Ap  (r)  r2  dr  (5o) 

di  Jo  uo 
=  -Q  (r) 


Httr" m SS&C ££c£iS£2££iLjsa  x  tl-%  a  f 


<*r~  T~iV~~i~i-i  t‘\  T't 


«wssc*?j<» 


where  Q(r)  is  the  total  charge  displaced  inside  a  sphere  of  radius 


r.  Thus  we  find 


*£  =  -.-a £si 

dr  t,<~ 


(5-h) 


Hence,  the  electric  field  in  the  Kanzaki  lattice  (neglecting 
the  P1  terns  in  eqn.  (  5-1  ))is  given  by 


£  (r)  =  -  V$ 

displaced  charge 


=  r  Q(r) 

r  r2 


(5.5) 


Adding  on  the  electric  Held  of  the  ion,  of  resultant  charge 


Z  equal  to  the  valency,  we  find 


e  (r )  =  r  ["q(i*)  -  ‘l  \  • 


(5-6) 


For  free  electrons,  we  know  asymptotically  that  we  have  tho 
Eriedel  oscillations  represented  by 


Apuo  ~  Acos  2kf  r 


(5-7) 


where  k^  is  the  Fermi  wave  number.  Wo  can  therefore  write  for 
the  asymptotic  form  of  Q(r)  -  2 


r* 

ihe  introduction  of  a  phase  shift,  necessary  for  strong  scattering, 
is  easily  effected. 
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Q(r)  -  Z 


-/. 


OO 


APuo  to 


r°° 

■4*  A/  cos  t  dt 

J2kfr  t 

4*  A  Ci  (2kf  r) 


(5.8) 


where  GL(x)  is  the  usual  cosine  integral,  as  shown.  Hence  we  nay 
write 


'I 


£  (*)  ~  r  4*  A  Ci(2kf  r) .  (5<9) 

r4~ 

An  estinate  of  A  is  readily  available  freer  the  work  of  Stott, 
Baranovsky  and  March  (1970).  Sane  results  for  Q(r)  in  Cu  are 
shown  in  Fig  3.  Hence,  if  we  assume  this  electric  field, 
evaluated  at  the  atonic  positions  in  the  Kanzaki  lattice,  acts 
on  the  resultant  ionic  charge  Ze  at  each  site,  we  have  a  first 
estinate  of  the  Kanzaki  forces  as  Ze  £(r). 


We  want  to  emphasize  that  this  is  only  a  very  crude  example . 

It  will  be  necessary  to  estinate  the  contribution  due  to  p  ^  in 
eqn.  (  5.1  )  as  well  as  to  calculate  ApUQ  beyond  a  f reo-ele ctron 

model.  In  both  calculations,  we  need  the  perfect  lattice 
solutions.  Work  i\  in  progress  to  evaluate  the  displaced  charge 

j 

^P-qo  in  II  for  a  mono-vacancy  in  the  unrelaxed  lattice,  from  ; 

the  theory  of  Beeby  (1967),  but  results  are  not,  as  yet,  available.  ^ 

J. 

t 

However,  Harris  (1968)  has  already  successfully  applied  Beeby* s  j 

> 

\ 

nethod  to  tfc9  impurity  problem.  Since  the  Beeby  theory  is  based  \ 

t 

on  the  KKR  nethod,  the  necessary  response  functions  are  already  J 


1 

.  tj 


at  afea 


•1 

-26- 

contained  in  the  KKR  solutions  for  the  perfect  lattice.  Thus 
we  have  here  a  basis  for  a  rather  precise  evaluation  of  the 
displaced  charge  in  the  relaxed  defect  lattice. 

Once  the  electric  field  is  known,  the  Kanzaki  forces  can  be 
found  and  Kanzaki *s  original  jaethod  used  to  obtain  a  new  set  of 
displacements. 


1 


'Vt  |k*V 


*  v®5*-.-  S  S-^St  v^ 
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6.  Kanzaki  forces  and  energies  of  mono-  and  divacancy 


As  a  further  illustration  of  the  present  approach*  let  us 
try  to  get  the  Kan  said,  forces  associated  with  a  Hartree  treatment 
of  a  vacancy  in  a  metal.  Then  we  can  write  for  the  total  energy 
of  the  metal 


U  = 


+ 


❖ 


is) 


+  * 


as  p  (r)  p(e) 


(6.1) 


The  first  term  in  eqn.  (6.1  )  is  the  kinetic  energy,  the  second 

is  the  classical  electron-electron  interaction  energy,  tho  third  is 
the  interaction  energy  of  the  conduction  electrons  with  the  ion 
cores,  while  $  represents  the  core-core  interaction  energies. 


As  the  present  writers  have  argued,  the  energy  can  ha  written 
when  the  metal  is  deformed,  in  terms  of  the  components  of  the 
displacement  Held  as 

U  =  U°  +  *  L  A^'aP  ^  ^*0  (6.2) 

Z  Z'  ap 

and  if,  at  the  same  time,  the  lattice  is  subject  to  external _ trees, 

then  the  following  term 

K  “  ~  Yj  P^a  (6-3) 

Za 

must  be  added  to  eqn.  (  6.2  )  F  are  the  Kanzaki  force 

components. 


. . 


s5*Bi 


write 


’■"*v*f  ^v££,*>  »*^*4r?fv^t-*  ey%~£iz' 


- *■  '--.s' 


If  we  assvrae  the  existence  of  a  pseud  oaten,  then  we  can 

write 

-£  T?/f^-Cg~  ^pseudo  a  tom  +  Jdr  A(r)  JjV(r)-  t(r)J 

Act  i-  ^ 

-/dr  pv(r)  *(r) 


(6.ii) 


where 


*(r)  =  )  #(*  -  £)  •  16.5) 

It  will  clearly  be  of  Interest  to  see  whether  the  above 
approximation  to  the  Kanzaki  forces  gives  results  similar  to  those 
of  the  previous  section,  based  on  direct  calculation  of  the 
electric  field  at  the  ions. 


The  interest  iff  the  above  calculation  resides  in  its  possible 
extension  to  estimate  directly  the  effects  of  relaxation  on  the 
divacancy  binding  energy.  Unfortunately,  it  is  then  necessary 
not  only  to  know  how  the  displaced  charge  round  the  two-centre 

ft. 

problem  Is  related  to  that  round  the  one-centre  problem  ,  but 
also  it  is  necessary  to  relate  the  positions  of  the  relaxed 
atoms  round  a  monovacancy  to  thoss  when  the  second  vacancy  is 
brought  up  to  the  near-neighbour  distance. 


We  want  to  conclude,  because  of  the  difficulty  of  this  problem, 
by  pointing  out  a  possible  approach  related  again  to  calculating 
the  electric  field  E(r)  in  the  final  equilibrium  configuration. 

ft. 

Such  a  calculation  for  an  impurity-vacancy  complex  in  a  metal  was 
carried  out  by  Alfred  and  March  (1957).  For  a  diva  can  cy,  it  is 
practicable  to  find  the  displaced  charge  from  the  model  of  Seeger  and 
Bross  (1956). 


>>W.a  ...at 


We  conclude  that,  unless  the  use  of  pair  potentials  can  be 
justified,  which  is  unlikely  round  a  divacancy  in  a  polyvalent  metal, 
it  is  go?  ng  to  be  of  considerable  interest  to  map  out  the  electric 

field  E  (r)  (zero  at  r  ■  £+  u,  )  in  the  final  state,  and  find  the 

—  ~  —  — 

associated  energy  stored  in  the  field.  We  suspect  that  (  o  7  ) 

ought  eventually  to  be  replaced  by  difference  equations  in  an 
exact  formulation. 
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We  turn  finally  to  discuss  how  the  energy  In  the  Kanzaki 
lattice  might  be  calculated  within  the  present  framework. 


Let  Tg  and  UK  be  the  kinetic  and  potential  energies  in  the 
Kanzaki  lattice,  for  a  given  separation  R  between  the  two  vacancies. 
Ihen  we  can  write 


a  Tv  +  'Jv  =  -  \ 


L^  +  '- ie>- 


(6.9) 


where  the  P  are  the  Kanzaki  forces.  Naturally,  those  forces, 
and  also  the  displacements  Ug  depend  ei  the  separation  R. 


In  one-body  potential  theory  as  used  in  this  paper,  where 
an  attempt  is  made  to  incorporate  some  of  the  many-body  effects 
into  the  potential,  it  seems  best  to  relate  the  tf»u*l  energy  to 
the  potential  energy  since  at  least  the  classical  part  of  this 
is  calculable  from  the  electron  density  plus  the  potential  of  the 
nuclear  framework.  Ihus  we  can  write  the  above  ecuation  in  the  form 


VK>  +  VR)  *  -4  (-  +  a*)*  ^ 


(6.10) 


Evidently,  fran  the  electron  density  and  the  electric  field  e(r) 
wo  can  calculate  E^(R)  and3J^(  co  },  the  difference  giving  us  a 
part  of  the  relaxation  energy.  Naturally  E^(  no  )  involves 
crucially  the  displacements  around  a  single  vacancy,  while  E^R) 
involves  those  round  a  divacancy. 
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If  X-ray  scattering  experiments  could  be  refined  to  yield 
not  only  the  electron  density  pQ(r)  in  the  perfect  crystal,  but 
the  density  around  an  impurity  or  an  imperfection,  band  theory 
would  allow  the  construction  of  the  crystal  potential  V  (r*y  and 
the  defect  potential  Vj(r).  Until  such  time  as  that  becomes 
possible,  we  must  use  the  best  available  methods,  based  on 
refinements  of  Slater  p  ^  exchange  to  construct  these  potentials 
so  as  to  incorporate  at  least  some  part  of  the  electron-electron 
interactions. 


Knowing  starting  displacements  u^,  ,  (r)  can  be  estimated 

seml-empirically  from  eqn.  (  2.1  )  at  the  Fermi  level.  Further¬ 

more  its  energy  dependence  is  accessible  through  the  approximate 
result  (  U.lh  )  in  terms  of  the  local  denrity  of  states  in  the 
perfect  crystal.  The  charge  displaced  in  the  Kanzaki  lattice  by 
the  defect  potential  V^Qr)  can  then  be  estimated  from  that 
displaced  by  Vj  inserted  in  tha  unrelaxed  lattice,  corrected  by 
eqn.  (  li.l6  ).  This  displaced  charge  can  then  be  used  to  find 
the  Kanzaki  forces  and  hence  to  calculate  a  new  set  of  displacements. 


-a 

f 

% 


Finally,  it  is  emphasised  that,  within  the  present  framework, 
calculations  of  defect  energies  should  be  attempts a  from  the 
classical  potential  energy  terras,  plus,  if  necessary,  estimates 
of  the  sxchango  energy  from  the  electron  density.  The  kinetic 
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energy  should  not  be  calculated,  however,  from  the  density  matrices 
discussed  in  this  report,  for  these  have  incorrect  off-diagonal 
elements,  though  the  diagonal  elements  agree  v~  ih  the  correct 
many-body  density  matrix. 


Part  B 


7 .  Introduction 


In  earlier  reports  (see  also  Rousseau,  Stoddart  and  March, 
1970;  1971)  wo  developed  a  density  matrix  theory  of  electron 
states  in  disordered  systems  and  we  applied  it  specifically  to 
calculate  the  electronic  density  of  states  in  liquid  Be.  Our 
conclusion  was  that  the  dip  in  the  density  of  states  in  crystalline 
Be,  resulting  fran  the  band  overlap  which  explains  its  metallic 
character,  is  only  partially  > filled  in*  on  mnlting.  Vfe  argued 
that,  to  describe  this  effect,  we  needed  a  strong  scattering 
theory,  and  that  was  developed  in  our  earlier  wcrk. 


To  bring  such  a  theory  into  direct  contact  with  e-'porlment 
is  at  present  difficult.  Although*  in  principle.  photoemission 
experiments  on  liquid  Be  ought  to  reflect  the  density  of  states, 
it  is  true  at  present  that  the  most  cirect  contact  with  experiment 
is  through  the  electrical  transport  coefficients.  !<fe  therefore 
develop  here  a  theory  of  the  electrical  resistivity.  The 
motivation  for  the  present  approach  is  provided  by  the  theory  of 
inverse  transport  coefficients,  pioneered  by  Edwards  (1965). 

This  theory  is  therefore  discussed  first,  in  its  general  form, 
and  subsequently  single-centre  approximations  to  the  density  matrix 
are  made,  in  order  to  make  progress  in  the  calculation  of  the 
resistivity.  Finally,  a  method  of  calculating  single-centre 
density  matrices  fran  radial  wave  functions  is  presented,  humeri. cal 
calculations  based  on  this  approach  are  in  progress. 


% 


f^.4^^}  *  '  -X«Vv»  ^$5 


8.  An  exact  formula  for  the  electrical  resistivity 


The  derivation  which  follows  closely  resembles  Greenwood* s 


(19.58 )  treatment  which  led  to  an  exact  expression  for  the 


conductivity.  The  final  fom  for  the  resistivity  bears  a  strong 


resemblance  to  an  expression  obtained  by  Edwards  (1965).  We  show 
that  Edwards  *  result  needs  modifying,  however,  in  the  -denominator 


term. 


Consider  a  Hamiltonian 


2  2  2 
H(t)  =  1_  )  8_  +ePt\  +(%  b_  )  +Ai  ,3^1 

2m  L  \  i  0x  /  \  i  3y/  V  i  32/  J 


(8.1) 


+  #(x+ct,y,z) 


where  the  vector  potential  (-eFt,0,0)  corresponds  to  an  electric 
field  (F,0,0)  and  the  potentials  of  the  system,  giving  rise  to 


the  total  potential  $  as 


(E)  - 1 


(8.2) 


are  all  moving  with  velocity  (-c,0,0).  If  c  is  small,  then  we 


nay  write 


and  so 


$(x+ct)  z  $(x)  +  ct  j3  $(x) 

*3x 


H(t)  =  H_(t)  +  ct  3$ 
°  3x 


(8.U) 


yqgassffwggg 
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where  Hq  is  got  by  putting  c  -  0  in  (8.1)  ,  and  is  in  fact  equal 
to  the  Hamiltonian  used  by  Greenwood  (1958).  He  proves  that 

HQ(t)  “  exp(-ieFtx/&)Ho(0)exp(ieFtx/&)  (8.5) 

and  so  from 

H(t)  «  exp(-iePtx$i)H(0)exp(ieFtxyii)  +  ct  $  (8.6) 

where  jS>*  =...  3$/3x,an8  H(0)  is  the  Hamiltonian  at  t  *  0, 
which  is  the  Hamiltonian  of  the  unperturbed  system.  Obviously, 

H  will  have  an  eigenvalue  equation  of  the  term 

H(t)Y(t)  =  E(t)iT(t)  (8.?) 

and  if  one  takes 

?n(t)  =  exp(-d.ePtx/4i)  (8.8) 

tdiere  ^  »  ^(0),  can  easily  be  seen  that  the  eigenvalues  in  (8.7) 
are  equal  to 

En(t)  -  En  ♦  t*^  (6.9) 

where 

$nn  =  &)  ♦»&)  (8.10) 

But  ^*nn  as  defined  in  (8.103  is  just  the  expectation  value 
of  the  force  on  an  electron  in  the  state  n  in  the  absence  of 
any  electric  field,  and  hence  is  zero. 

The  equation  of  motion  for  the  density  matrix  p  is 

dfi.  =  Cgp] 
at 


(8.11) 


3S"gr=rf 


■5hv^-<;  •■/tT~/?'*~*-Xf?i-,« 


Substituting  this  into  the  equation  of  motion  tor  the  roatrix  elements 


of  p  ,  i.e. 


J  P^  =  (n||aW  +  (|S  |pSb)*NpJ-  ,  (8a2) 


and  using  the  standard  result 


(m  |  dn)  =  -(m  |  dH  1  n) 

«  _av,„... 


(8.13) 


we  arrive  at  the  result 


P  -  f 
Hnn  n 

where  fR  ■  f(E^),  the  Pbmd.  function  for  energy  En,  and 


(8. Hi) 


%  ■  +  c**nJ  (f§)E  6<W  (8-«) 


where  we  have  written  P  -  f  +  g,  and  worked  to  first  order  in  F 
or  c  throughout.  Also,  » 


*  %  [dr  ♦"  (r)  3  ^_vr) 

im  J  n  Bx  m 


(8.16) 


i.e.  is  the  velocity  matrix  element. 


The  Hamiltonian  which  wo  have  used  (equation  (8.1)) 
corresponds  to  a  physical  situation  where  the  coordinate  exes 
are  translated  with  uniform  velocity  c  down  the  electric  field. 
If  we  make  c  equal  to  the  drift  velocity,  then  we  are  in  a 
position  to  calculate  the  force  on  an  electron,  which  of  course 


V~>> 


rSptyiJ ?'$'^5^>^>'>»*- '»•«! /r~'v "&S  ??  ^a»"  " 
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must  vanish,  since  a  steady  state  has  been  established  hr 
and  Hence 

-/  j  ^nm  ®  nn  +  e^)  pnyt  =  0  (8.17) 

nm  4  ii 

for  c  equal  to  the  drift  velocity.  Neglecting  oscillating  terms, 
we  got  the  resistivity  R  to  be 


R  =  -2^1 


2  2 
ftp  e  run 

v? 


E  ••-  (i)»  • 

ran  n 


«(B-Sn) 


(8.18) 


_  -I 

fl,+2^  Vv  (W\  S(E-E)| 

l  tyT  L>  ^  \cLE/E  m  j 


Eqn  (8.18)  can  be  written  in  terms  of  the  energy  derivative  of  the  Dirac 
matrix  p  *  where 


p’(rr'E)  =  ^  ♦nV)'Jfn(£’)6(E-EI1) 


and  the  result  is 


R  -  N/(1+») 


(8.19) 


(6.20) 


where 


N  = 


-2t& 

*pI 


le2/“aEf/^^I^)I>)  *>2i> 

a  ^ 


I  p’fe^g  b)J 


and 


D  =  Z4x~ 
fiimp0 


f  dS  df  f  dr^o  9£(£2)  p' <2,^)3  p’C^E) 

Jo  313  J  dx2  0X1 

(8.22) 


where  II  is  the  total  volume,  and  p  the  mean  electron  density. 

w 

The  tern  N  is,  apart  from  multiplying  factors,  the  same  as  Edwards 
(1965)  finds,  but  the  denominator  term  D  is  different,  in  that  it 
contains  P ’ ( r2  E),  whereas  Edwards*  result  replaces  this  by 
G+  (£^£2  ®)»  It  is  not  cccrroct  to  say  (as  Edwards  has  done)  that 
G+  contains  terms  involving  [  E-E^I  which  are  related  to  the 
virial.  Rirther,  it  is  easy  to  see  that  (8.22Xs  real,  whereas 
Edwards  *  denominator  would  give  rise  co  a  complex  resistivity. 


Hie  exact  repression  for  the  resistivity  is  given  by  (8.20)  and 
is  the  basic  result  in  this  section.  We  proceed  now  with  an 
analysis  of  this  formula,  to  obtain  a  foxm  which  is  amenable  to 
numerical  evaluation.  From  this  point  we  work  in  atomic  units 
e  »  %  «*  m  “  1. 


Representation  for  the  Dirac  matrix  p  ’ 


To  enabl  che  ensemble  average  to  be  performed  in  equation 
it  is  useful  to  consider  a  representation  for  the  Dirac  matrix  in  the 


p'(£+x,£»S)  =  P’(x,E)[[  ^1+f(£"ii»  £»E)) 


(9-1) 


Substitution  of  (8.1)  into  the  numerator  N  of  (8.20) 
gives  after  sane  re -arrangement  of  factors. 


<N>  = 


%i /®  «Ea  (p’0<^-iia»E))‘ 


(10.1) 


k\{\  (' + 


t£rV  E>))  a=9=o 


where . 


fc<e(^i>E2*E)  =  (1+f(^£2E))(1+f<£c£1B5) 


3#(ij)  +  9  -1 


r  °»»V£ 

1°“^ 


(10.2) 


The  major  difficulty  in  the  evaluation  of  the  ensemble 
average  in  (10.1)  is  lack  of  knowledge  of  the  n-boqy  ionic 
correlation  functions,  only  the  two-bo^'v  correlation  function  being 
reasonably  well  known. 


We  may  obtain  immediately  the  result  for  a  random  system 
from  (10.1)  as  in  this  case  the  correlation  functions  - 
are  of  course  known  oxactly.  Wo  obtain  (sec  .Rausseau,  Stoddart 
and  March,  19?0) 

2 

^:>R  =  JdEft(r,)/rdv1dr2  (pHrrv2,7.)'j  (10.3) 


where 


average  <  N>  ^  takes  the  same  general  form  as  the  random 
average  result  (10.3)  with  Fr  replaced  by  a  function  of 
■which  contains  information  about  the  ionic  correlation  through 
the  two-body  correlation  function,  if  we  use  the  approximation  s 
proposed  by  Rousseau,  Stoddart  and  1-brch.  Even  if  the  exact 
correlation  functions  were  used  to  perform  the  average,  the 


-lo¬ 
re  suit  for  <  N  >  would  retain  the  sane  structure  as  (10.3) 
with  F^  replaced  by  a  function  involving  these  exact  correlation 
functions.  However,  even  the  crudest  approximation  to  these 
functions  seems  to  lead  to  an  intractable  result  for  <  N  > 
and  we  will  consider  here  only  the  random  average  result  (10.3). 
We  propose  later  that  an  approximate  account  of  structure  can 
be  taken  simply  through  the  liquid  structure  factor  S(g) 
in  a  manner  identical  to  that  in  the  usual  pseudopotential 
formulation  (Ziman,  1961). 


11.  The  Weak  Scattering  limit 


To  0(0  ;  the  denominator  in  (8.20)  does  not 

contribute  and  we  obtain  from  (10.3)  and  (10.  U) 

2 


e 


s=~2n  /“*’(*>  jdr 1  tog  42)  <WE> 

(n.i) 

PR(2)  •*?/«*  (VrT  (V^<£e -?->)•  (31.2) 


where 


If  to  define  the  Fourier  transform  0  (f)  by 


e(r)  =  1  >-  ~ 


i  r  •-  -  /  i 

J 


(11.3) 


and  use 


p»  (jl,-e2.e)  =1. 


j  1  2  6(E-:(k2/2» 


(n.U) 


i 


& 


then  we  obtain. 

(r,-£2,E) 


and  through 


1  E 


{£i-^2)  |?(B)|2 


(11.5) 


o 


R  = 


l^Pe 


I  dq 


43  }?  (4)1  2 


(11.6) 


This  result  is  the  well  known  Ziinan  formula  for  the  resistivity 
if  we  regard  $  (<j)  as  the  Fourier  transform  of  the  pseudopoteniial 
describing  the  interaction  between  an  electron  and  a  single  screened 
ion.  It  may  be  shown  that  an  equivalent  result  can  be  obtained 
from  the  numerator  N  in  the  liquid  metal  case,  the  only  change 
being  the  replacement  of  the  potential  part  of  (11.6)  with 
S(q)  $  %  (q)  }2  where  S(q)  is  the  structure  factor. 

12.  Definition  of  Pseudopotential 


Jp-, 


The  single  ion  potential  cannot  be  regarded  as  a  perturbation 
and  we  require  an  evaluation  of  the  complete  term  (’10.3)  .  We 
work  with  the  random  average  result  which  can  bo  put  into  an 

rs» 

interesting  form  if  we  define  F(k,E)  by. 


m 


¥ 


fr(^-£5»s)  =_1 


ife. 


F_  (k,E)  •  (12.1) 


E3' 


fe 


I 


wk 


R  = 


*/. 


2k, 


3  p  /d4  qJ  A- (a). 

12*V  o  * 

r  e 


(12.3) 


To  second  order  in  0  we  have  of  course, 

A(2)  (4)  =f#(q)|2 


(12.ii) 


but  in  the  general  case  (lO.Ii)can  be  used  to  define  AR(q).  We 
obtain  an  equivalent  result  for  the  liquid  metal,  the  contribution 
of  the  numerator  taking  exactly  the  same  form  as  (12.3)  .  However 
in  this  case  the  appropriate  transform  A^(q)  depends  on  both  the 
local  ionic  potential  0  and  the  ionic  correlation  functions. 


We  have  therefore  the  interesting  result  that  if  the  num¬ 
erator  N  is  the  dominant  contribution  then  the  resistance  takes 
exactly  the  same  form  as  the  weak  scattering  result  (11.6) 

We  may  use  A^(q)  therefore  to  define  a  local  ’pseudopotential 1 
V(r)  through 


Aj  (4)  =  |  V(q)  | 


2 


(12.5) 
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in  the  usual  sense  as  that  potential  which  gives  the  correct 
scattering  phase  shifts  in  Born  approximation.  In  the  liquid 
metal  cans  it  is  obvious  that  this  * pseudopotential 1  will  depend 
on  the  structure,  although  in  a  first  approximation  one  may  take 
the  potential  V(r)  from  the  random  average  result  and  include 
structure  simply  by  use  of  A^(q)S(q)  in  place  of  A^.  A  structure 
‘pseudopotential'  arises  in  the  formulation  of  Rubio  (1969)  and 
has  been  calculated  by  Ashcroft  and  Schaich  (1970).  However 
Rubio » u  result  for  R  is  correct  only  to  lowest  order  in  (kf  &)”1 
where  ^  is  the  mean  free  path  of  an  electron  at  the  lersd  surface 
The  -theory  presonted  here  is  free  from  such  limitations. 


13. The  Denominator  Term 


The  ensemble  average  of  D  given  by  (8.22)  can  be  performed 
in  a  similar  way  to  that  of  N.  However  we  use  *ere  the 
additional  representation, 


V  (£y£g>S)  3po  (&,.-££*£)!  j  -t^rg-t^E) 


ox. 


8xi 


a 


(13-1) 


(13.2) 


The  final  result  for  the  random  average  is 
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<  E  >R  -  -  g^lfi.  J jj,  f^Ej/dr.drg/a^  pjUj^-Z.Sg-Z.E) 

pen  '  "'  3x, 

II 

x  9p£  teg~£grj^  (13.3) 

3x. 

t 

where  terms  of  higher  order  in  p  have  been  omitted.  To  obtain 
the  correet  ensemble  average  for  R  we  must  of  course  average  both 
N  and  D  together.  However  a  sufficiently  good  approximation  will 
probably  be  given  by 


<  R‘  >=  <  N  >  * 

1+<  D  > 


(13.U) 


The  interesting  possibility  then  arises  that  in  certain  circum¬ 
stances  (if  the  scattering  is  sufficiently  strong)  wo  could  get 
1  +  <  D  >  _»  0,  representing  a  transition  to  an  insulating 

state.  However  this  point  requires  further  investigation,  end 
we  will  focus  attention  here  on  the  contribution  of  N  to  the 
resistance. 


U4.  Proposed  Computatio>^l  Scheme 

We  propose  a  computational  scheme  based  upon  ihe  result  (12.3) 
for  the  resistance,  but  including  the  effect  of  structure  through 
the  introduction  of  the  liquid  structure  factor  S(q).  Thus  we 
have  approximately, 

*  An  improved  approximation  is  <K>  ®  <N>  -  <ND>  which 

TT^TS  > 

reduces  to  (13. U)  if  <ND  >  -  <N>  <D  >  • 

s 


(lii.l) 
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2kf 

R  =  p  f  dq.  q3  A^q)  S(q) 

lS?Pe2  ■'o 

in  whi^h  therefore  the  random  ensemble  average  of  N  is  used  to 
define  a  ‘pseudo-potential*  according  to  (12.5)  Ihe  structure 
enters  through  S(q)  in  a  manner  identical  to  that  in  the  standard 
pseudcpotential  form  for  R. 

A  further  simplification  is  assumed  in  the  formula 
for  the  resistance  in  the  first  calculations.  We  obtain  AR(q) 
from  the  Fourier  transform  of  the  first  term  only  in  (10. lr) 
neglecting  the  effect  of  the  exponential  tern.  Thns  using  (10.5) 
we  take 

V£T£2’E)  =  P  ]*&  <K£r->)  (V-r20(^2"^)  (lU,2) 


vhere  AR(q)  is  determined  by  (12.1)  and  (12.2)  ,  Equation  (lb-2) 
is  just  the  term  of  0( p )  in  an  expansion  of  FR  in  powers  of  the 
ionic  density.  What  we  are  essentially  assuming  is  that  the 
‘pseudopot antial*,  which  is  defined  in  terms  of  AR  by  (12.5) 
does  not  have  any  explicit  dependence'  oh  the  ionic  density. 

We  turn  finally  to  discuss  the  single-centre  scattering,  defined 
through 


Hie  potential  is  assumed  to  be  spherically  syinmetrlc. 


15.1  Tha  t-matrix  Its  integral  equation  is 

tCXflO  =  v(x)6 (x-£)  +Jdz  v(x)  Gq(x  -z)t(z,£)  (X£#1) 

v(x)  is  the  single  centre  potential  and 
+  . 

Go  (£“£')  =  “(4x)  exp[+ik  \  r-r’{  ][1  r-r *J  3"*1  (15*2) 
i  x 

•where  k  *  E  ,  E>  0,  k  “  i(-E)  ,  E<  0.  Defining  a  2-sided  Fourier 
transform  of  t  thus : 

*(£’£)  =  j&x  d£  exp[i£.x-i  a.£]t(x,£)  £^3) 

and 

t^»y)  s  (2*)~6j&R  d£  exp[-i£.x  +  i£.jr]t(£,q) 
and  appropriately  transforming  (l£.l)  the  xeeult  is 


Writing 


*(£,£)  =  ▼(£-£)  +  (27c)"3Jdk  v(£-k)G0(k)t(k,£)  (15.5) 

l  * 

«**.<Vp)  W w  (15.6) 


and  similarly  for  /.cj  v(g-^)  the  individual  terms  in  the  expansion 
of  (15. 5)  are 

■tCpq.)  =  v(pq)  +  (2x)  3  J  k20k  v(pk)GQ(k)t(kq)  (15.7) 

+ 

The  suffix  6  has  been  dropped  for  convenience.  0o”  (k)  can 
be  written 


^(k)  =  [E-k2]-1,  -ix6(E-k2) 

V/  T 


(15.8) 


(16.1) 


16.  Solution  of  the  t  matrix  for  continuum  states  (E>0) 


Separate  t  into  its  real  and  imaginary  parts: 

t(pq.)  =  R(pq)  +  il(pq) 

Substituting  (16.1)  and  (15.8)  into  (l£.7)  ve  get 


R+(p<l)  =  v(pq)*(4*)  2v(p^)vEl+(vl3Q)+/dkQ(pkE)rf'(kq)  (16.2) 

00  o 

I+(pq)  =  -(4*)~2v(p/Ey£  R+(vEq)+J  a.k  G(pkE)l+(kq)  (16.3) 

-1  ° 

where  C(pkE)  “  (2x  )~^  v(pk)k2  [E-k2  ]  .  Equations  (16.2)  and  (16.3) 

are  of  standard  Fjedholm  type,  and  are  solved  by  determination 
of  the  inverse  kernel  to  G(pkE),  which  satisfies  the  following 
integral  equation: 

00 

F(pkE)  =  G(pfcE)  +  J  dk'  G(pk"E)F(k'kE)  (16.U) 

o 

In  terms  of  F,  the  solutions  of  (16.2)  and  (16.3)  are 

K+(pq)  =  w(pq)  +  oj(p/EVE  I  (vEq)  (4x)  “2  (16.5) 


I+(pq.)  =  -w(pyE)  Ve  R(VSq)  (4*)' 


(16.6) 


where 


oi(pq)  =  v(pq)  +  /  F(pkE)  v(kq) 


(16.7) 


Equations  (16.5)  and  (16.6;  can  be  separated  by  putting  p  °\/E  in  (16.6) 
substituting  into(l6.5)  and  then  putting  p  -  V  E  in  (16.5)  • 

This  gives  R(  -/Eq)  solely  in  terns  of  (0  ,  and  this  expression 
for  R(\/Eq)  can  be  substituted  back  into(l6.6).  The  results  are? 
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£ 


I 


m 


1 


m 


R+(pq)  =  w(p<l5“(4x)”^  E  a)(p/B)o>(v5VS)w(Vteq) 

[l  +  (4*)~^  E^v/EVt;)2]-1 

(16.8) 

*  (po)  =  **(4x)  2  V&  w(py/E)a)(vIlq) 

[t  +  (4*)“^  ^(Vivlil)2]”1 

(16.9) 

It  is  easily  seen  that  R*(pq)  »  R+(pq),  I**(pq)  ■  -I+(pq). 
17.  The  solution  for  p*(r  r'E) 


Sie  double  BsUrlor  transform  of  v(x)  6  (k-y)  according  to  eqn  (15.3) 


is 


v(£-.a)  =  j  dx  exp[i(£-£).x]  v(x) 


(17.1) 


Expanding  both  sides  of  equation  (l?.l)  in  spherical  harmonics 
it  is  easy  to  see  chat; 


CO 

V(pq)  =  (4*)2  J  x2  dx  ^(>.)3(px)  j(qx)  (17.2) 


o 


so  v(pq)  "  v(qp).  Expanding  to(pq)  in  a  perturbation  series 
by  means  of  equations  (16.70  and  (16 .It)  it  can  bo  soon  that  every  tern 
is  symmetric  with  respect  to  interchange  of  p  and'q.  •  Hence, 
because  of  equation  (17.2). 

w(pq)  =  to(qp)  (17.3) 

Hence  also,  frcm  (16.0)  and  (16.9)  R(pq)  “  R(qp),  I(pq)  "  I(qp)* 


The  relation  between  the  t  matrix  and  •  the  Green  function 


is: 


°0b£)  =  %  (*»£)  +  j^SP^  %  ; )%  '— 

J  (17 • U ) 


whore  G_is  given  by  equation  (15.3).  Double  Fburler  transform  of 


yields: 


G(pq)  =  (2*)3  6(p-q)p“2  G.(p)  +  G0(p)t(£#a)  Gtf(fl) 


(17-5) 


Expanding (17 . 5 Hr.  spherical  harmonics  using  equation 


the  terms  are: 


G(pq)  =  (2*)3  6(p-q)p”2  <%  (p)  +  Gd(p)t(pq)GQ(q) 


Using  the  standard  result  that 


(17.6) 


p* (pqE)  =  iCar)"1  [G*(pc|E)  ~G~(pqE)]  (17.7) 


and  substituting  (16.8)  nod  (16.9)  into  (17.6),  tfao  result  is 


*p'(pgE)  s  (2*)3  6(p-q)p"’2  7c6(E-p2) 

-  m2w(pv)3M3co(\/Eq)6(E-p2)6(E-q2)[l+(^)*”UE“(v5>/3)^^) 
+4w<pq)-(U«)^w(p/]3)^>^)tt(^qXl+(^)^^  )2]  3 

x[6(E-p2)(E-q2)“1  +  6(E-q2)(E-p2)“1] 

+  (Ux) ~2VE<»(vSEl)  o)(n/EJ q)  [  1  +  (4*) oWm)  2  ]-1  [  (B-p^ (E-q2 )  l”1 

(17.8) 

Pn^rn  the  expansion  of  equation^. Unapplied  to  pO  in  spherical 


harmonics  it  is  seen  that: 


(2k) **(4h) ~2p '  (xr 'E)=  yp^dpq^ dqj (pr )  j (qr • f ) p ’ (pqE) 


t  2,__2. 


(17.9) 


}* - ‘Z *4$!S&j3 ■»V>^ v^'^-- =5?3 «?«  *5  • 


19.Friedel  sun  rule  and  idenrpotency  of  Dirac  matrix 


The  Dirac  matrix  is  obtained  from  p’  by  an  integration 


over  energy 


E 

p(r,r’E)  =[  dE*  p’(r,r’ ,S)  (19.1) 


It  is  idempotent,  i.e. 


P(£.£',E)  =  Jdx  p(r,x,E)p(x,r‘,  E)  (19.2) 


and  obeys  the  closure  relation 

p  (r  ,r.’ ,  oo)  =  &(r-r '  ) 

In  equation  (19.1)  only  the  continuum  states  are  considered. 
The  bound  state  contribution,  ?iz. 

P-U (£»£.’)  =  ^^*(£)  0i(£l) 


(19.3) 


(19. h) 


must  be  added  in  before  equation  (19.3)  is  true.  Here  i  runs 
over  all  the  bound  states  of  the  system.  If  the  scattering 
potential  is  to  be  completely  screened  at  large  distances,  then 
the  displaced  charge  must  equal  the  net  charge  producing  the 
potential,  i.e. 


jd Z  tp(  rr 


E.n)  -  Prt(r,r,E.)]  =  Z 


(19.5) 


where  E^  is  the  energy  of  the  highest  occupied  state,  and  the 
region  of  integration  is  large.  Expanding  equations  (19.1)  to  (19-5) 
in  spherical  harmonics  it  is  seen  that 


■fiafe 


$6- 


E 

p(rr’E)  =  J  &Bk  p*  (r  r’E) 


(19.6) 


p(rr‘E)  =  J” x2dx  p(rxE)p(xr‘E) 


(19.7) 


n 

f>(rr*ao)  s  r  6(r  -  r‘) 


019.8) 


2.Ux  V  (21+1)  /“  R  r2dr[p(rrEf)  -  pQ(rrE  )]  =  Z  (19*9) 

4*  J 


The  factor  of  2  in  equation  (19. 9)accounts  for  spin  degeneracy. 
The  radial  vave  function  for  continuum  states  satisfies  the 
di  '  otial  equation 


u.  r 
dr2 


R(kr)  -jj^k2  ~  v(r)  J  r  R(kr)  =  O  (19*- 


vhere  E  =*  k  and  rR(kr)  °  0  at  r  »  0.  Multiplying  ^19.10).  by 
cos  r)  (k),  setting  vn  a  similar  equation  for  rR(k*r),  subtracting, 
integrating  by  parts  gives 


D’co8T)(k,)rR(k,r)l  &_  [cos-q(k)rR(kr)3  -  [cosT)(k)rR(kr) 

dr 

d.  [co3T\(kt )  rR(k*r)”j 

dr  J 

=  (kr2-k2)  J  R  r2d r  cosT^kJcosi^k’ )R(kr)R(k’r) 

0  (19.11) 


£ 


Si 


fe 


gf 


p 


£ 


I 
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u 
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I 
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Bie  asymptotic  form  of  R(kr)  is 


R(kr)  =  j(kr)  -  tanrj(k)  n(kr) 
Hence,  substituting  j(  p)-+  .p"1  Cos  [p*4(^+l)*3 , 
n(p)  =  P  sin[p-4-(-6+l  )tc]  and  multiplying 


(19.12) 


by  cosrj(k)  gives  the  result 


coa  T](k)R(kr)  =  (kr)-1'  Pinfkr  4  rj(k)  -  &%)  r  -  oo  (19.13) 


Substitution  into(19.H)  gives 


/V 


r  dr  coa'q^cosT^k*  )R(kr)R(k’r)  =  i 

•  01.0..  t  /l. 


2kk T  (k+k  ’) 
viiich  is  the  same  thing  as 


2kk’(k-k‘) 
sin[  (k-k*  }R+Ti(k)-Tj(k' )]  (19.1h) 

sin[(k+k’)R  +  r,(k)  +t|(k’)  -  C%] 


£ 


r  dr  cos'^(k)cosTl(k,  )R(kr)R(k'r) 


r 


2kkT  L 


6(k-k*  JcostT^k)^^  )1 


-C(k+k’)cos[-n(k)+T)(k,)-^TC] I  +  coa(k<»k,,)R  sin[T}(k)-T](k’) 

J  2kk’(Jc-k’) 


-  cosCk^k*  )R  sin[,n(lc)+TiO,'T  )-&ft] 
2kk5(k+k’) 


(19.15) 
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uslng  the  representation  6(x)  =  dim  sin  gx/’xx,-  . The 

g-*oo 

result  for  free  electrons  Obtains  by  setting  "H  m  0* 


Hence 


I  r2dr[p '(rr'k)-p^(rr 'k)] 


=■  _JL  8inrn(k)"^i(kM1  +  an  oscillating  term* 

2k2  k’**k  k-k»  (19.16) 

But  since  ‘n(k*)  =  t)(k)  +  (k'-k)  drj/dk  and  neglecting  the 


.2  k’-k 


oscillatory  term: 


^  r2dv[p  *  (rr *k)-p  (rr‘k)]  =  dTT/|k  =  i  drj  (19.17) 

A  A1_  2  -  dS  w  ' ' 


2k  y/E 


Hence  using  (18.9)  and  (19.9)  we  obtain  the  result 


*  2^(2d+!)  Tl(Ef)  =  Z 


(19.18) 


which  is  just  the  Friedel  sum  rule.  To  prove  that  the  solution  (17.15) 


produces  an  idempotent  Dirac  matrix,  using  (19.6)  and  (19.7)*e  obtaiu 


p(rrfE)  =1,  /  2k1?.  dk,2k,,ik"coc^l(k,  )cosTl(k,,)R(rk, ) 

~2  Jo 

*  R(r’k") 

X  />  cosn(k')oosT,(lc")R(x!c,)P.(xk"  )  (19.19) 


£'»■  ■*<?..  *^X*»  »-  <>  l^,.-  ,VW 


which  is  to  be  proved,  if  the  solution  is  idempotent.  Hie 
last  factor  in  equation  (19. 19 )is  given  by  equation  (19.15).  Notice 
first  that  the  terra  in  6  (k»+kn)  doss  not  contribute,  since  k-*' 
and  k"  are  always  greater  than  0.  Secondly,  the  oscillatory 
terms  also  give  zero  contribution  as  R  -*»  ,  since  the  integrands 

are  smooth  functions  of  k*  and  k".  The  only  part  where  care 
is  required  is  the  vanishing  denominator  in 


coafk-k^R  si 
2kk’. 


when  k’ 


k-k' 


But  as  (19.6)  cud  (191?)  chow,  this  Gives'  rise  to  no  difficulty.  Hence 
00 

[  x2dxcos■^(k,)co8TJ(k,,)R(xl:,)R(xk,,)  =  %  6(k-k')  (19.20) 

Jo  2k2 

and  so  the  right  hand  side  of  equation(19.19)  becomes 

r  k  p 

1  /  2k*  dk’  oosTiCk*  )cosT](k’ )  R(rk*  )R(x',k') 

*  Jo 

which  is  jusi  p(rr'E),  by  equations  (17 .15) and  (19*6).  Adding 
in  the  bound  state  contribution  to  P  gives  no  difficulty,  as 
this  part  is  idempotent  in  itself,  and  integrals  of  products  of  $ 
and  R  vanish,  as  these  functions  are  orthogonal. 


20.  Proof  of  closure  relation 


Expand  an  arbitrary  function  ^(r  r* )  iu  the  function 
P  ’(rr’E)  and  the  bound  state  wave  functions  0^(r)  (we 
write  P^(r)  for  the  ^th  component  of  ^(r),  for  simplicity). 


Then 


Ccrnporing  (20. U)  and  (20.5  )ve  see  that  we  must  have 


*  yy  _ 

f  2k  akcosir}(k)R(kr)cosT\(k)R(kr* ) 
K  JO 

=  aCr-r’l  -  V 

1  1 


(20.6) 


Hence  the  soxution  given  in  equation  (18.9)  together  with  the  bound 
state  contribution,  satisfies  the  closure  relation,  as  well  as 
idempotency  and  the  Priodel  ana  rule. 


I 
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21.  Nummary 

Part  A  has  reported  a  treatment  of  relaxations  round  point 
defects  in  metals,  using  electron  theory.  Preliminary  results 
are  reported  for  Kg  in  la  and  for  a  vacancy  in  Cu,  and  further 
work  is  in  progress  on  applications  of  this  theory. 

In  Part  B,  our  earlier  work  on  electron  states  is  extended 
to  deal  with  the  electrical  resistivity  of  liquid  metals  where 
there  is  strong  electron  scattering.  The  theory  should  be 
applicable  to  liquid  metals  like  calcium  and  barium,  which  have 
relatively  high  resistivity  and  numerical  calculations  on  calcium 
are  in  progress,  using  the  method  outlined  immediately  above  for 
treating  the  single-centre  scattering  from  a  calcium  ion. 
Generalizations  to  liquid  metal  alloys  are  also  under  investigation. 


Ve  want  to  summarize  here  some  of  the  basic  properties  of 
the  functions  F,  R  and  jP  which,  in  principle,  can  be  found  from 

the  theory  of  the  perfedt  crystal.  To  obtain  F,  as  stressed 

i 

above,  is  a  one-body  prdblem,  once  the  potential  V^(r)  is  known. 
On  the  other  hand,  exact  determination  of  R  and  P  is  not  possible 
at  prese.  't  as  these  quantities  depend  on  the  exchange  and 
correlation  energy  of  an  inhomogeneous  electron  gas. 


Proof  of  relation 


)  between  response 


function  F  and  one-body  Green  function  and  densii 


matrix 


As  remarked  above,  the  relation  (  3.1*  )  was  given  by 

Stoddart,  March  and  Stott  (1969).  An  elementary  proof  of  this 
will  be  sketched  below. 


Bin.  (  3.7  )  gives  us  the  change  in  the  Bloch  density 

matrix  C  (r  r*  B  )  due  to  a  change  in  potential  AV,  a  change 

o  —  — 

we  denote  by  AC(r  8  )  on  the  diagonal. 


Then  wo  have  immediately  from  oqn.  (  3.7  )  that 


-■^v'  -r\>-x~:  *\"'  - 


AC(r  ^) 


o 


<^•1  CoCiir’P  -  P1)  C0(r‘r  31). 

(Al.l) 


We  nw  cctapare  this  result  -with  the  density  change  ^P  (r  E) 
caused  ty  the  sane  change  in  potential.  This,  from  eqn.  (1*-H  ) 
jd-th  Vd(r)  set  equal  to  zero  is 

r 

Ap(r  E)  =  Jar1  AV(r')  F(r  r*  B).  (A1.2) 


But  ve  have  the  Laplace  transform  relation  between  C  and  p 
(provided  we  add  a  positive  potential  energy  to  bring  all 
energies  £  0) 

AC(r0)  =  pj  Ap  (r  B)  e"PE  aE.  (M-3) 

J  o 

Integrating  this  equation  by  parts,  we  find 


d_  Ap(rE)dE 
dm 


=  jar1  AV(r’)  j°°  e"pE  0F(r  _r'  E)  dS 
0  3E 

where  the  last  line  follows  by  differentiating  eqn.  (  A1.2 
respect  to  E. 


(Al.l*) 


)  with 


It  is  clear  then  fren  eqns  (  Al.l  )  and  (Al.l*  )  that 

3F/3E  can  be  related  to  the  integral  over  g  of  the  product 

1 

of  C0’s  displayed  explicitly  in  eqn.  (  Al.l  ),  Substituting 


aBSasaa 


<■  4*  ‘»4&+WM  ^*H(V^v 
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the  explicit  fora  (  3.1  )  for  CQ,  the  3^  integration  is  readily 

accomplished.  Similarly,  the  inverse  Laplace  transform  required 
to  obtain  oF/  3  E  from  equ.  (  Al.li  /  is  easily  completed.  Hence, 
using  the  explicit  definitions  of  CQ  and  pQ,  eqn.  C  3.1i  )  follows. 

Integral  equation  for  R 


Hie  determination  of  R  from  first  principles  presents  more 
difficulty,  though  a  send-empirical  procedure  can  be  devised  to 
obtain  a  useful  sts-ting  approximation  (see  Appendix  3). 


The  argument  sketched  below  follows  the  discussion  of  JK 
and  is  given  for  completeness  in  that  the  potential  Av(r)  required 
to  generate  the  Kanzaki  lattice  depends  on  P  through  eqn.  (3.2  )• 

Hie  essential  point  is  that  we  can  write  the  change  in  the  one- 
body  potential  as  we  move  the  ions  from  j&  tod  +u^  as 


AV  =  AV 


elect  ostatic  +/u(r  r’)p1(r’)  dp’.  (41.5) 


But  p^  is  related  to  R  (r)  and  we  also  have 

<1 


AV  =  y  k,  R  (r*  -  A)ap-  .  Y,  [ 
e  Jr  -  r* )  t 

+  ^  u^.  Ju(r  r’ )  R(r*  -6)  dr1 


-  Ze 

5*  TF=4 


(ai.6) 


1  v<"vv5''- r  •"  'r:- ' -'  ~^^ZT' 


•*'t****'Af+'t*  •v'-S'-f’*  •'V’JsV 


T~Cr  r^*-  V iV‘/V'"w’'--‘  J-'1 
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Ihis  can  now  bo  expanded  to  0(u  ^  )  to  yield  the  form  (■  3.2  ) 

where 


F.{^'  )dr 1 

is  “  ft 


+ jTJ (r  vj  )  R  (r  *  )dr 1 
r  (A1.7) 


After  sene  manipulation,  the  basic  integral  equation  (3.3  ) 
of  JM  then  follows. 


i 


l 

$ 


\ 

3 

l 

i 

4 


fe-rT.  wi~^ 
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1 


1 
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Appendix  2 


To  avoid  the  troublesome  problem  of  solving  eqn.  (  ^.2  ) 
by  numerical  iteration,  we  shall  now  argue  as  follows.  With 
the  defect  potential  Vj(r)  introduced  into  the  perfect  lattice, 
let  us  suppose  that  by  the  methods  of  Stoddart,  March  and  Stott 


(1969)  we  have  generated  a  density  matrix  Cu,  there  u  stands  for 


the  unrelaxed  defect  lattice.  Then  we  define 


AC  -  C  -  C 
uo  u  o 


and 


LCfK  *  Cf  “  V 


(A2.1) 

(A2.2) 


It  seems  reasonable  to  assume  that  a  useful  starting  point 


for  the  calculation  of  Ac_.  would  be  A  C  ,  which,  at  least, 

Ziv  uO 


includes  the  full  effect  of  the  defect  potential  Vj(r),  though 
not,  at  first,  the  relaxations. 


It  is  now  a  simple  matter  to  show,  from  the  integral  form 
of  the  BOLoch  equation,  that 

iCfK  =  -/  P  aPl./dr"[(4CfK<£  P  -  Pl>  *  °fK(£  E’P-Pt’) 


Vd(r")  CK(r,,r’  (A2.3) 


and 


^i3AL -_JL:  -w>— 


St^r^ycays^iaq!^  frtV*!-"  v^ts:  rk**.  4<  ;, 


*fVi-,A^  >4  ***»v.>t  _*t  j^sc's* a. f  a^ST 
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4Cuo  “  -j^,JA^'  [  (4f!uo(s.  E"P-Pt)  »  Guo(r  r"0-p.,)) 


Vd(£u)  C0(r“  r*  p.,). 


(A2.h) 


Subtracting  these,  vie  have  for  the  error  dC  =  AC^  -  AC^q 
the  result 
0 


5C  =  -J  dp^dr1'  6C  +  ACU0  •.•  CKJvd(r")CK(rM  r*^) 


'/o  ^l/^"  t°w>  +  °.]^>  C°(-  ~  &1} 


(A2.5) 


Now  we  have  also 


°K-Co*  °K0 


(A2.6) 


whore  Aj&>  is  0(u  ^  ).  Thus,  to  0(u£  )  we  find 

P  r. 


w  =  -j  ^1jdrH  ac  +  ACKoJvd(r")  Cjr"  r'^.) 


-/  ^tf^T  [A  Cuo  +  Co]  ^Ko^'Pl^ 


(A2.7) 


If  we  further  assures,  as  starting  point,  that  6C  <AC. 


Ko 


then  we  have  as  the  fir3t  approximation  to  correct  the  charge 
•displaced*  by  Vd£r.|)  in  the  unrelaxed  crystal 


60  ~  -2 f  V  dp1  Jar”  ACKq  Vd(r,‘)  CQ(r"  r*  p^ 


(A2.8) 


-  /  B  a  fit  jar-  iOuo  Va  <r»)  acKo  (r»  r-  pt). 
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Here  then,  is  an  expression  to  correct  the  displaced  charge 
for  the  effects  of  relaxation,  in  terns  of  the  properties  of  the 
unrelaxod  crystal  containing  the  defect  potential  V^r). 

To  get  a  self-consistent  theory,  it  is  clear  that  eqn.  ( I^.io  ) 
must  be  used  to  refine  the  displaced  charge  in  the  unrelaxed  crystal. 
If  we  work  with  the  Green  function,  the  (3  integration  in  ( ii.10  ) 

is  immediately  removed. 

It  may  well  be  possible  to  calculate  local  changes  around 
the  vacancy  by  making  the  additional  approximation  in  eqn.  ( U.10  ) 
that  varies  slowly  in  space.  Then  for  the  diagonal  element 
6G(r  r  p)  we  find 

6C(r  r  p)  « 

-  va(p) J  00  4B1  far"  [4CK0  (r  r'  (3,) 

^2Cq  (r"  r  P  -  P1 )  (&• 9 ) 

+  AC3uo(r"  r  P  -  P^J 

In  the  main  report  we  give  &  rather  cruder  fora  of  eqn.  (  AH. 9  ). 

Clearly,  however,  eqn.  (  A2.7  )  provides  a  perfectly  proper 

starting  point  for  refining  the  displaced  charge  calculated 


in  the  unrelaxed  lattice 


■&** Wfc!?’  ,<-3LW<f‘JM  .rtiw*  ■jt'ti  t.-tfJt-pFi _ 
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Appendix  3 

Preliminary  estimate  of  p^ (r),  the  charge  'displaced1 
in  the  Kanzaki  lattice,  for  vacancy  in  Cu 

Since  the  evaluation  of  the  displaced  charge  in  the 
presence  of  relaxation  is,  via  eqn.  (  li.ll  ),  intimately 
related  to  the  charge  displaced  in  Kansaki  lattice,  namely 
P^(r),  we  have  made  scrae  preliminary  estimates  of  this 
quantity  for  a  vacancy  in  Cu. 

Evidently,  to  calculate  p^(r)  fran  eqn.  (  2.1  ), 

we  must  have  information  on  the  displacements  and  the 
response  function  R  (r). 

Tewordt  (1958)  has  dealt  with  the  problem  of  the 
relaxations  round  a  vacancy  in  Cu  by  matching  local  relaxations 
around  the  vacancy,  obtained  using  Bo to -Mayer  pair  potentials, 
on  to  the  long-range  relaxations  which  can  be  properly  treated 
by  replacing  the  discrete  crystal  by  a  continuum.  This  approx¬ 
imation  beccnes  valid  outside  a  sufficiently  large  region 
enclosing  the  defect  or  defect  complex.  For  then  the 
displacements  u  ^  become  small  and  vary  slowly  from  atom  to 
at  cm. 


We  have  used  these  displacements  u  ^  as  given  by 
Tewordt  (3.958),  even  though  they  will  eventually  need 
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re  fining  to  take  account  of  the  detailed  nature  of  the 
displaced  charge  round  the  vacant  site  as  discussed  in  this 
report. 

It  only  then  remains  to  set  up  a  suitable  approximation 
to  R(r)  in  eqn.  (  2.1  ).  We  already  have  the  Fourier 

components  at  -the  reciprocal  lattice  vectors  K  through 
eqn.  (  2.3  ).  Bitterman,  Chipman  and  de  Marco  ((1961) 

have  studied  the  X-ray  scattering  from  Cu  and  we  therefore  have 
approximate  estimates  available  of  the  pK*s.  Following  JM, 
we  have  made  the  assumption  that  we  can  use  the  fom  (  2.3  ) 
for  all  k,  smoothly  interpolating  between  the  pK*s  at  the 
reciprocal  lattice  vectors  to  obtain  p  ^ .  Obviously  this 
procedure  is  not  utlque  and  eventually  the  integral  equation 
of  JM  (3ee  Appendix  1)  must  be  solved  to  find  R(r). 

Hie  above  assumption,  as  emphasized  by  JM,  is  equivalent 
to  the  assumption  of  pair  forces,  and  we  can  then  write  R 
as  the  gradient  of  a  scalar  density  a (r),  namely 

R  (r)  =  7o(r)  (A3.1) 

where,  from  eqn.  (  ?.2  )  it  then  follows  that 

P0(r)  =  ^  o(r  -  £).  (A3.2) 

£ 

Provided  o(r)  has  the  correct  Fourier  components  at  the 
reciprocal  lattice  vectors,  which  is  of  course  ensured  by  our 
procedure  to  within  the  accuracy  of  the  experimental  measurements. 
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then  oqn.  (  A3. 2  )  will  repr*  -duco  the  exact  gi-ound-otate 

charge  density  of  the  crystal. 

The  form  of  o  used  to  fit  the  X-ray  results  of  Batterman 
at  al  is  given  (cf.  Jones,  March  and  Tucker,  1965}  in  Figs  Al  -  A3. 


placed  charge  in  Cu  metal 


As  an  example,  we  shall  compute  the  displaced  charge  in 
Cu  metal  using  Tevordt '  s  displacements  and  making  the  prelim¬ 
inary  approximation  R  °  Vo  .  We  fit  o  the  ‘pseudoatan1 
charge  density  to  the  X-ray  scattering  experimentsas  shown. 
in  Figs  Al  and  A?. 


Using  the  Hartree-Fock  atom  density  as  starting  point, 
the  correction  Ao  (k)  required  to  fit  the  X-ray  scattering 
at  the  Bragg  reflections  is  shown  explicitly  in  Fig  Al. 

This  is  not  unique;  we  have  simply  drawn  a  smooth  curve 
through  the  measurements,  and,  of  course,  eventually  the 
integral  equation  of  Jones  and  March  will  have  to  be  used  to 
find  R(r).  Hie  Fourier  transform  Ac  (r)‘  ie.  shown  in  Fig  A2. 


* 


Vo(r  -  £) 


(A3* 3) 
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for  various  directions  in  the  crystal  and  the  results  are 
shown  in  Figp.  A3~Af> .  The  very  marked  anisotropy  oTidant 

there  is  a  consequence  of  the  fact  that  as  we  go  out  along 
the<  111>  direction  we  do  not  encounter  an  atom  over  the 
range  shown.  This  then  is  a  first  estimate  of  the  charge 
displaced  in  the  Kanzaki  lattice,  strained  to  account  for 
the  displacements  round  a  vacancy  in  Cu. 


Appendix  U 


Relation  between  response  function  approach  and  elasticity 
theory  in  long  wavelength  Unit 


The  customary  elasticity  approach  gives  us  the  result 
that,  in  response  to  a  localized  body  force,  which  in  k  space, 
in  the  long  wavelength  limit  gives  us  the  Fourier  transform 
of  the  force  F  in  the  form 

F  (k)  “  ik  X  (Ah.l) 

where  X  measures  the  strength  of  the  body-force,  the 

* 

•  displacements  (r)  take  the  form 


(r )  =  c  r 


(AU.2) 


The  multiplying  constant  in  eqn.  (  Ah.2  )  is  determined  in 
Tewordt 1 s  method  by  matching  to  the  ‘discrete*  displacements 
in  the  immediate  neighbourhood  of  tne  defect. 


It  is  evident  in  Kanzaki‘s  method  that  F  and  u  are 
related  by  the  dynamical  matrix  Dap  (k),  which  in  turn, 
as  <JM  discuss,  is  related  to  the  response  function  R(r). 

It  then  becomes  clear  that  elasticity  theory  gives  us  rather 
direct  information  on  the  response  function  R(k)  at  small 
k,  and,  in  particular  relates  the  small  k  behaviour  rather 
directly  to  the  elastic  constants.  This  will  be  discussed 
quantitatively  elsewhere  by  Ciaesson,  Jones  and  March  (to  be 
published)  in  connection  with  many-body  forces  in  lattice 
dynamics. 
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We  want  only  to  make  one  other  point  relating  to 
elasticity  theory  here.  Corresponding  to  the  displacements 
(  Al;«2  ),  jt  is  well  known  that  there  is  a  volure  change 
per  defect  g'.ven  by 

Aft  =  Uxor  (Ab.3) 

idiere  Y  =  3(1  -  °)/(l  +  o),  a  being  Poisson*s  ratio. 
W®s.anticipatQ  that  eventually  it  nay  be  necessary  to  check 
that  the  defect  potential  adopted  obeys  some  condition  equivalent 
to  a  Frledel  sum  rule  into  which  the  volume  change  (  Ah.3  ) 
is  incorporated,  as  we  indicated  in  toe  Introduction.  Our 
microscopic  theory  is  not  yet  sufficiently  well  developed  to 
give  a  precise  understanding  of  this  point,  which  is  worth 
further  study. 
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Captions  for  Figures 


Pig.  1.  Local  density  of  states  0(r  E)  for  Mg  in  Id  (solid 
curve)  for  E  *»  |  E^.  Dashed  curves  show  local  density  of  states 
in  pure  Mg  and  pure  Id  metal. 

Pig.  2.  As  in  Pig.  1  but  for  E  a  E^private  cozanunication  from  J. 
Inglasfield). 

Fig.  3*  Total  displaced  charge  Q(r)  inside  sphere  of  radius 
r  in  case  of  vacancies  in  at  (2  ■  -1),  Mg(Z  =  -2),  A -A  (2  •  -3) 
and  Pb  (Z  -  -U). 

Ibis  is  needed  to  calculate  electric  field  in  Kanzaki 
lattice.  Results  3hown  are  from  the  work  of  Stott,  Baranovsky 
and  March  (1970;  see  also  March  and  Stoddart  (1968;  p.?5l) ) . 

rt  o  ^ 

Fig.  A.?..  Ao  versus  (h  +  k  +  ^  )*  for  Tu,  used  to  fit  X-ray 
results  of  Batterman  et  al.  Form  is  gi/en  by 

Ao(K)  =  +  A2K7  e  d2K 

A1  «  3.2h,  a1  =  l.li,JU  »  0.CO83,  =  1.3?. 

Pig.  A. 2.  A o  (r)  corresponding  to  Ao  (K)  in  Pig.  A.l. 

Fig.  A. 3.  Charge  density  pQ  in  Ca,  together  with  displaced  charge  P 

in  Kanzaki  lattice  when  strained  by  a  vacancy.  Densities  shown  along 
< 100>  direction. 

Pig.  A.lt.  Same  as  Pig.  A.3.  hut  along  <  110 >  direction. 

Pig.  A.?.  Same  as  Pig.  A.3.  and  A. It.  but  along  <  111  >  direction. 


Q  Along  ^lOO^*  Direction 


Fig  A4 


b  Along^lO?’-  Direction 


Fig  A5 


c  A!ong<1!1>  Direction 


